The Spin Glass Phase in the Four-State, Three-Dimensional Potts Model 
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We perform numerical simulations, including parallel tempering, on the Potts glass model with 
binary random quenched couplings using the JANUS application-oriented computer. We find and 
characterize a glassy transition, estimating the location of the transition and the value of the critical 
exponents. We show that there is no ferromagnetic transition in a large temperature range around 
the glassy critical temperature. We also compare our results with those obtained recently on the 
"random permutation" Potts glass. 
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I. INTRODUCTION 

Potts models are among the building blocks of sta- 
tistical mechanics, and their disordered versions (Potts 
glasses) are commonly used to describe a large class 
of anisotropic oricntational glasses.^ For example, if a 
crystal of molecular nitrogen is disordered by includ- 
ing some percentage of argon, the resulting compound, 
Ari^x{^2)x, is a disordered quadrupolar glass. S The 
Potts glass is one of the models of choice to describe 
materials of this type. 

The four state {p = 4) pure Potts model in two di- 
mensions {D = 2) describes the adsorption of N2 on Kr 
in graphite layersi^ In D = 3 it describes the behavior 
of FCC antiferromagnetic materials (NdSb, NdAs and 
CeAs for example) where the magnetic field points in 
the <1,1,1> direction^. The site diluted version of the 
Potts model can describe, for example, the adsorption of 
hydrogen on the (1,1,1) plane of nickel which has been 
previously disordered by inserting oxygen atoms. ^ In this 
paper we will study the four state glassy Potts model, 
in which quenched random disorder induces frustration. 
This model presents at least three interesting theoretical 
problems that are still unsolved. 

The first of these is the nature of the spin glass phase 
transition; one needs to reliably compute the critical tem- 
perature and the critical exponents in order to character- 
ize the universality class of the model. 

The second issue is how the qualitative features of 
phase diagram, including spin glass and ferromagnetic 
phases, evolve when going from the mean field models to 
realistic, finite dimensional models. For example, previ- 



ous work^i^ has shown that, at low temperatures, the 
standard Potts glass (in which the Potts coupling be- 
tween two spins can have positive or negative sign) de- 
velops spontaneous ferromagnetic ordering, which can af- 
fect, or even prevent, a spin glass phase transition. Fur- 
thermore, in mean field theory the value of this ferro- 
magnetic transition temperature, Tfm, varies with the 
number p of Potts states as Tfm = (p — 2)/2, which gives 
Tfm = 1 for the four state model. Mean field analysis 
also shows that for p < 4 the spin glass transition temper- 
ature (where replica symmetry gets broken) is Trsb — 1: 
if p = 4 the two transition temperatures coincide. Hence 
an interesting open problem is to check whether or not 
this result, valid for D = 00, also holds when the di- 
mensionality is finite. Also relevant here is that, in the 
mean field picture, the p = 4 glassy model is "marginal" , 
since for p < 4 the transition is continuous whereas for 
p > 4 the order parameter q(x) is discontinuous (even if, 
as usual in spin glasses, there is no latent heat)^^ In a se- 
ries of interesting papers Binder, Brangian and Kob^ (see 
also the recent work by Katzgraber, Lee and YoungiS) 
study the ten state glassy Potts model, and find that for 
such high number of states the mean field and the finite 
dimensional cases are very different. Here we investigate 
whether the same is true for p = 4. 

The third relevant issue is again related to universal- 
ity. In order to avoid a possible contamination of the spin 
glass transition point by the effects of the ferromagnetic 
phase, Marinari, Mossa and Parisi^^ have introduced a 
new class of glassy Potts models, the "random permuta- 
tion" Potts glass, where a gauge symmetry protects the 
model against a ferromagnetic transition. This approach 
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has the advantage of being closer to reahty since, in real 
quadrupolar glasses, ferromagnetism plays no role. One 
of these models has been thoroughly studied recently by 
some of the authors of the present work,^^ and its be- 
havior found to be consistent with a Kostcrlitz-Thoulcss 
phase transition. One of its signatures is that, below 
the critical point, data for the correlation length divided 
by lattice size, ^/i, for different sizes merge into a sin- 
gle curve. However, given the precision of the numerical 
data and allowing for corrections to finite-size scaling, 
one cannot exclude a value of the lower critical behavior 
near and slightly below D — 3. 

A further motivation for this study is to investigate 
how the behavior of the Potts glass changes with p. For 
p = 3, a Potts glass transition occurs^S, with critical ex- 
ponents v ~ 1.2 and rj ~ 0.02, while, for p = 10, Ref. [l^ 
finds no phase transition in agreement with Ref. [§]. In 
addition, we note the Ising spin glass model, which cor- 
responds to p = 2, has I' ~ 2.5 and 77 ~ —0.4 (see 
Ref. (UHIl). It would therefore be very interesting to 
get a consistent picture of how the nature of the Potts 
glass transition evolves with the number of states p. 

In an attempt to give reliable answers to these ques- 
tions, we have performed extensive numerical simulations 
using one unit of the JANUS dedicated computer (which 
has a total of 16 units). — We have been able to thermal- 
ize the p — A Potts glass model on an L = 16 cubic lattice 
down to the low temperature phase: this gives us infor- 
mation on far larger lattice sizes than had been possible 
before. 

The outline of the paper is the following. In section 
2 we introduce our model and physical observables. In 
section 3 we describe the numerical methods that we have 
used in the simulations. In section 4 we describe our tests 
of thermalization and our approach to data analysis, and 
we analyze our findings, both for the overlap and for the 
magnetization. The main results are that we have been 
able to characterize the spin glass transition and that no 
onset of ferromagnetic order has been found at and below 
the spin glass transition point. We discuss these findings 
in section 5. 



II. MODEL AND OBSERVABLES 

In the p-state Potts model, each site z of a three dimen- 
sional cubic lattice of linear size L with periodic bound- 
ary conditions has a scalar spin Si which takes one of the 
values 1, 2, . . . ,p. The Hamiltonian of the standard Potts 
glass model is 

H = -Y, .h, Ss^,s, , (1) 

where the sum runs over all nearest neighboring sites. 
Two neighboring sites i and j interact with energy —Jij 
when their spin states Si and Sj have the same value, 
and otherwise their energy is zero. The couplings Jij 



are independent quenched random variables taken from a 
bimodal distribution {Jij = ±1) with zero average. From 
now on we will focus on the four-state {p — 4) case. 

It is possible to represent the state of site i by a (p — 1)- 
dimcnsional vector. Si, equal to one of the p unit vectors 
Sa pointing to the corners of a hyper-tetrahedron in (p — 
l)-dimensional space. These vectors satisfy the relations 

S.-S, = P-^. (2) 

p- 1 

Equation ([2]) defines the simplex representation, that we 
will use to describe the observables measured in the sim- 
ulations. 

In order to investigate the possible presence of (spuri- 
ous) ferromagnetic effects we have carefully checked the 
value of the magnetization, looking for possible signs of 
spontaneous ferromagnetic ordering. In the simplex rep- 
resentation we define the vector magnetization as 

1 ^ 

i=l 

where N = is the number of spins. 

The existence of a possible transition to a ferromag- 
netic phase has been also analyzed by studying the mag- 
netic susceptibility 

XM = Nj}^, (4) 

where ((• • •)) stands for the thermal average and (• • •) 
denotes the disorder average. 

To study the glass transition we define the spin-glass 
order parameter as a tensorial overlap between two repli- 
cas (i.e. two copies of the system defined with identical 
couplings whose spin values evolve independently). The 
standard definition of its Fourier transform with wave 
vector k if^- 

'i"'ik) = j^T.S^'^"^^"^'^'""'''' (5) 

i 

where S^^^^ is the /i-th component of the spin of the first 

replica in the simplex representation, and S^'^'''^ is the 
I'-th component of the spin in the second replica. 

The momentum-space, spin-glass susceptibility is de- 
fined by 

X,{k) = NY^W^^im . (6) 

We also define the correlation length ^ in terms of the 
Fourier transform^^ in Eq. ([6]) as 

I f^-lY' (7) 

where km is the minimum wave vector allowed within 
the lattice. Periodic boundary conditions imply that this 
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vector is fern = (27r/L, 0, 0) or one of the two other related 
vectors in which the components are permuted. The def- 
inition in Eq. ([7]) arises naturaUy on a finite lattice. 

We will base a large part of our analysis on the dimcn- 
sionless correlation length ^/L, i.e. on the correlation 
length divided by the lattice size. This quantity does not 
depend on L (asymptotically for large L) at the tran- 
sition temperature, which allows us to obtain a precise 
estimate of from the value of T at which data for dif- 
ferent lattice sizes cross.— 



III. NUMERICAL METHODS 

We have simulated three dimensional cubic lattices 
with linear sizes L = 4, 6, 8 and 16. Because spin-glass 
simulations have very long relaxation times, we used the 
parallel tempering (PT) algorithniii to speed up the dy- 
namical process that brings the system to thermal equi- 
librium and eventually explores it. Physical quantities 
are only measured after the system has been brought to 
equilibrium. 

The dynamics is comprised of single-spin updates and 
temperature swaps. The single-spin updates are carried 
out with a sequential heat bath (HB) algorithm. We 
define a Monte Carlo sweep (MCS) as N sequential trial 
updates of the HB algorithm (i.e. every spin undergoes 
a trial update once). 

The PT algorithm (applied to a given realization of the 
quenched disorder, that we will call a sample) is based 
on simulating a number of copies of the system with dif- 
ferent values of the temperature but the same interac- 
tions. Exchanging the temperature of two copies with 
adjacent temperatures with a probability that respects 
the detailed balance condition is the crucial mechanism 
of PT. The result is that each copy of the system drifts 
in the whole allowed temperature range (that has been 
decided a priori) . When a copy is at a high temperature 
it equilibrates fast and so each time it descends to low 
temperature it is likely to be in a different valley in the 
energy landscape. 

The HB and the PT algorithms require high quality 
random numbers; we generate them with a 32-bit Parisi- 
Rapuano shift registerJ^ pseudo-random number genera- 
tor. 

Details about our numerical simulations are summa- 
rized in Table m The simulation of the smaller lattices, 
with L — A and 6, was performed on standard comput- 
ers. More powerful computational resources are needed 
to deal with the L = 8 and 16 systems, so we have studied 
them on a prototype board of the Janus^^ computer, an 
FPGA based computer optimized for a relatively small 
set of hard computational problems (among them, spin 
glass simulations). A performance comparison between 
an Intel(R) Core2Duo(TM) processor and one Janus pro- 
cessor (one FPGA) shows that the latter is about one 
thousand times faster— when simulating Potts models. 
JANUS has allowed us to thermalize a large number of 



samples for bigger sizes than would have been feasible 
on a standard computer. The computational effort be- 
hind our analysis amounts to approximately 6 years CPU 
time on a 2.4 GHz Intel(R) Core2Duo(TM) processors for 
L — 8 and thousands of CPU-years for L — 16. 

Data input and output is a critical issue for JANUS 
performance, so we had to carefully choose how often 
to read configuration data; in general, we end up tak- 
ing fewer measurements than in simulations on a tradi- 
tional PC. Having fewer (but less correlated) measure- 
ments does not affect the quality of our results. We read 
and analyze values of physical observables every 2 x 10^ 
MCS for both L = 8 and 16 (see Table |T] for details). 



L 


samples 


MCS 






iVHB 




4 


1000 


3.2 X 10^ 


[2.0,6.0] 


9 


5 


10^ 


6 


1000 


8 X 10^ 


[2.5.5.0] 


7 


5 


10^ 


8 


1000 


2 X 10** 


[2.7,4.2] 


16 


10 


2 X 10^ 


16 


1000 


8 X lO'' 


[1.7,4.1] 


32 


10 


2 X lO'^ 



TABLE L For each lattice size we show the number of disor- 
der samples that we have analyzed, the number of MCS per 
sample, the range of simulated inverse temperatures P — 1/T, 
the number of (uniformly distributed) /3 values used for PT, 
the number of MCS performed between two PT steps (A^hb), 
and the number of MCS between measurements {Nm)- 

On the larger lattices, we perform a PT step every 10 
MCS while on the smaller lattices this value is 5. In a 
standard computer the PT algorithm takes a negligible 
amount of time, compared to a whole MCS. However, 
in JANUS the clock cycles needed by one PT step are 
more than those needed for a MCS. For this reason we 
chose to increase the number of MCS between two PT 
steps. However, this number should not be too large, as 
we do not want to negatively affect the PT efficiency. A 
preliminary analysis has been carried over to test how the 
PT parameter would affect the simulation results, and we 
have selected a value that seems to be well optimized (see 
Table HJ. 

IV. RESULTS 
A. Thermalization Tests 

We start by briefly discussing the tests that we per- 
formed to check if our numerical data are really well 
thermalized. We use a standard test in which a given 
physical quantity is averaged (first over the thermal noise 
and then over the quenched disorder) over logarithmi- 
cally increasing time windows. Equilibrium is reached 
when successive values converge. We emphasize that it 
is crucial for time to be plotted on logarithmic scale. 

We are interested in the correlation length, defined in 
Eq. (O, which is plotted in Fig. [1] at the lowest simu- 
lated temperature (the hardest case for thermalization). 



4 



We see that the values of the correlation length reach a 
clear plateau for all sizes, strongly suggesting that our 
samples have reached thermal equilibrium. This anal- 
ysis also provides useful information about the number 
of sweeps that have to be discarded at the beginning of 
the Monte-Carlo history in order to use only equilibrated 
configurations. 



Applying Eq. (fTO|) to the operators dp^, and Xq yields 
respectively the critical exponents 1 -\- 1/v and 2 — 77^. 
Similarly, if we apply Eq. (fTO| to the magnetic suscepti- 
bility we obtain the exponent 2 — rjm. 



C. Overlap Critical Exponents 
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FIG. 1: A thermalization test. We show the behavior of the 
time dependent spin glass correlation length as a function of 
Monte Carlo time. We have averaged the correlation length 
using a logarithmic binning procedure. We show data for the 
lowest temperature simulated for each size. 




FIG. 2: The spin glass correlation length divided by L as a 
function of /? for L = 4, 6, 8 and 16. In the inset we magnify 
the crossing between the L — 8 and L — 16 curves. 



B. Finite Size Scaling Analysis; The Quotient 
Method 



To measure the critical exponents we used the quotient 
methodi ^^'^" In this approach one compares results for 
lattice sizes L and sL for integer s which here we take 
to be 2. Firstly, for a pair of lattice sizes L and sL, we 
find the point, (3 — /3cross, where the correlation length 
divided by system size is equal for the two sizes, i.e. 



sL 



or equivalently 



^{L, /3cross) 

L 



£,{sL, /Across) 
^{L, /Across) 



(8) 



(9) 



We then determine similar ratios for other observables. 
If an observable O diverges near the critical temperature 
as t~^° , where t is the reduced critical temperature, then 
we expect 

Qo{L,sL)^^^^^^l^^=s^o/^ + 0{L-n, (10) 

'-^[J^j Pcross ; 

where uj is the exponent describing the leading correc- 
tions to scaling. 



In Fig. [21 we plot the correlation length (defined in 
Eq. ([7])) divided by system size for different lattice sizes 
as a function of the temperature. According to Eq. ([H]), 
the data should cross if there is a transition. Indeed 
there are clear crossings, indicating a second order phase 
transition, though the values of (3 at which the crossings 
occur vary with system size. From Fig. [5] we determined 
the crossing values /3cross for the pairs of sizes (4, 8) and 
(8,16), see Table HIl By computing the spin glass sus- 
ceptibility and the derivative of the correlation length 
at these crossing points, we obtain estimates of the cor- 
responding effective critical exponents, rjq and v, from 
Eq. ([To]) and also show these results in Table [TTl Since 
we have only data at a discrete set of temperatures, we 
needed an accurate interpolating procedure to determine 
the crossing points and the values of other measurables at 
these points. We chose to fit all available data with a cu- 
bic spline. To test that our results are independent of the 
interpolation procedure we also implemented a linear in- 
terpolation aroimd the crossing point. We computed the 
crossing point and effective exponents with both proce- 
dures, and found agreement within the statistical preci- 
sion of our results. 

The two values of /Across shown in Table |lT] are rather 
different, suggesting large corrections to scaling, i.e. a 
small value for the correction exponent w, so we do not 
have enough information to reliably compute asymptotic 
critical exponents. Nonetheless, from Table |ll] we see 
that the trend of rjq with increasing size is very different 
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from what would be observed in the absence of a tran- 
sition for which r]g would equal 2. Hence, our numerical 
data strongly support the existence of a spin glass phase 
transition at finite temperature. 





/Across (-^1, L2) 




riq{Li,L2) 


r;m(Li, L2) 


(4,8) 


3.59(4) 


0.83(5) 


0.15(4) 


1.84(3) 


(8, 16) 


4.00(4) 


0.96(8) 


0.12(6) 


2.06(3) 



TABLE II: Results for the critical exponents using the quo- 
tient method. {Li, L2) are the two lattice sizes used and /3cross 
is the inverse temperature where the two curves of the dimen- 
sionless correlation length ^/L cross (see Fig. [2]). The values 
for V and r^q are extracted from measurements involving the 
overlap q, whereas rjm has been computed from the magne- 
tization. These results were obtained with the cubic spline 
interpolating procedure. 



D. The Magnetization in the Critical Region 




FIG. 3: We show the behavior of the magnetic suscepti- 
bility, xm, versus the inverse temperature. Notice that this 
susceptibility saturates in the critical region. In the inset we 
have plot ted the average of the modulus of the magnetization, 
< |m| >, against the temperature: this observable behaves as 
l/\/iV which clearly signals a paramagnetic behavior. 



As discussed in the introduction, the standard Potts 
glass studied here could undergo a ferromagnetic phase 
transition at low T. This second transition could bias 
our analysis of the spin glass phase by influencing the 
behavior even close to the glass transition temperature 
(a serious problem if the two temperature values are very 
close). It is therefore important to investigate whether 
there is a region with non-zero spontaneous magnetiza- 
tion close to the spin glass critical region. 

We have therefore computed, using the quotient 
method, the growth of the magnetic susceptibility, show- 
ing the results in the last column of Table [TTl The mag- 
netic susceptibility diverges with an exponent 2 — rjm , so 
rjm ~ 2 is a clear footprint for the absence of a magnetic 
phase transition. For the two largest lattices we find a 
value statistically compatible with 2. Hence we can safely 
discard the scenario where a ferromagnetic transition ap- 
pears at /Across- We are observing just a glass transition. 

In order to argue that there is no ferromagnetic tran- 
sition in the whole temperature range studied, we com- 
puted the magnetization and susceptibility throughout 
this range. In the paramagnetic phase, the magne- 
tization is random in sign so its modulus < |m| > is 
proportional to \/\/N, and the magnetic susceptibility 
Xm = iV< |mp > is independent of size. By contrast, in 
a ferromagnetic phase, < |m| > tends to a positive value 
at large N so xm diverges proportionally to N . 

In the main part of Fig. [3]we plot xm versus the inverse 
of the temperature. In the glass pseudo-critical region, 
(3 ^ 3.5 — 4, the two largest lattices give very similar 
results, so we recover the result 77^=2 quoted in Table [TTl 
Furthermore, at no temperature does the susceptibility 
increase strongly with size. Similarly, the magnetization, 
shown in the inset of Fig [31 decreases rapidly with size, 
which also indicates paramagnetic behavior. 



From Fig. [3lwe conclude that there is no ferromagnetic 
phase in the region /3 e [0, ~ 4.5]. 

V. CONCLUSIONS 

In this study we have numerically explored the equilib- 
rium behavior of a Potts glass with binary couplings on 
large lattices (L < 16 ). A prototype board (16 FPGA 
processors) of the Janusi^ optimized computer, using a 
parallel tempering algorithmic, has allowed us to do this 
for the first time. 

By computing the critical exponent associated with the 
magnetic susceptibility, and by analyzing the behavior of 
the magnetization in the critical region, we have shown 
that a paramagnetic-ferromagnetic phase transition does 
not occur. This result is different from mean field theory 
where, for general p, one sees both ferromagnetic and spin 
glass transitions at temperatures which become equal for 
p = A. 

We have found and characterized a spin-glass phase 
transition with critical exponents ~ 1, 77^ ~ 0.1 and 
hence (3q ~ 1/2. In order to extrapolate these values 
to the thermodynamic limit, larger lattice sizes need 
to be simulated (which we will try to accomplish in 
the near future). The critical exponents computed here 
are compatible with known values for other values of p. 
We note that the exponent v decreases with increasing 
number of states p, since v — 2.45(15) for p = 2 (see 
Ref. [3) and v = 1.18(5) for p = 3 (see Ref. [lO:]). 
The values presented in Table [TTl for p ~ A are consis- 
tent with this decrease which is expected to end when 
V = V = 2/D = (= 2/3 in Z? = 3), since a finite size 
scaling estimate implies that, for a disordered system, 
the transition is then first orderrSS Similarly, the expo- 
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nent t] grows with p since rj = —0.375(10) for p = 2 and 
r] = 0.02(2) for p = 3, while our estimates in Table HIl are 
larger. 

The hypothesis of a disordered first order phase tran- 
sition provides an upper bound rj — 1/2 (since the sus- 
ceptibility is expected to grow as L'^^^). In the mean 
field solution of Potts glass, second order phase transi- 
tions are found for small values {p < 4), but a first order 
transition is found^ for p > 4. An interesting problem 
for future study is whether the transition remains sec- 
ond order at large p for short-range spin glasses in three 
dimensions or whether, for a given value of p > 4, the 
second order transition disappears (in a tricritical point) 
to be replaced by a first order phase transition at larger 
pM 

We have studied the standard Potts glass which is ex- 
pected to be in the same universality class as the per- 
mutation Potts glass^^. However, the the present state 
of the art in numerical simulation does not enable us to 
confirm this. In the permutation Potts glass one does 
not observe a clear cut phase transition. Instead of a 
crossing point there is a smooth merging of the curves 
for different lattice size. This could indicate transient 
behavior, i.e. there is really a phase transition but it is 
only visible on larger lattices, or a Kosterlitz-Thouless 
like transition.™ For the standard Potts glass studied 



here, we find a finite transition as indicated by a crossing 
of the correlation length data in Fig. [21 However, we note 
that the crossing point shifts to larger /3, i.e. smaller T, 
at larger sizes. It is therefore possible that asymptotic 
critical behavior could be quite marginal, as is found in 
the permutation Potts glass. Consequently the standard 
and permutation Potts glass models may be in the same 
universality class, but larger sizes are needed to confirm 
this. 
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